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Radiative strength functions (RSFs) in 93_98 Mo have been extracted using the ( 3 He,a7) and 
( 3 He, 3 He'7) reactions. The RSFs are U-shaped as function of 7 energy with a minimum at around 



and 



E- 



3 MeV. The minimum values increase with neutron number due to the increase in the low- 



energy tail of the giant electric dipole resonance with nuclear deformation. The unexpected strong 
increase in strength below E 1 — 3 MeV, here called soft pole, is established for all 93-98 Mo isotopes. 
The soft pole is present at all initial excitation energies in the 5 — 8 MeV region. 

PACS number(s): 24.30. Gd, 24.10.Pa, 25.55.Hp, 27.60. +j 



I. INTRODUCTION 



The 7-decay of nuclei at high excitation energy tends 
to follow certain statistical rules. The dominating 7- 
transition driving factors are the number of accessible 
final states and the 7-ray transmission coefficient. The 
largest uncertainties are connected to the latter factor. 
In the description of this factor Blatt and Weisskopf [1] 
included an E^ L+1 dependency, where L is the angular 
momentum transfer in the transition. In their definition 
of the radiative strength function (RSF), this simple en- 
ergy dependence was divided out. With such a definition, 
the single-particle RSF (Weisskopf) estimates become in- 
dependent of 7-ray energy. Various concepts of RSFs and 
7 decay in the continuum are outlined in the reviews of 
Bartholomew et al. [2,3]. 

It has been well known that the RSF is not at all 
constant, but shows an additional E* dependency with 
x = 1 — 2 for 7 energies in the 4 — 8 MeV region. Axel 
[4] argued that this feature is due to the collective gi- 
ant electric dipole resonance (GEDR), which represents 
the essential mechanism for the 7-decay. However, the 
situation is more complex. Further studies [5-7] reveal 
fine structures in the RSF, which are commonly called 
pygmy resonances. This name does not refer to specific 
structures: the El pygmy resonance in the E y = 5 — 7 
MeV region of gold to lead nuclei could be due to neu- 
tron skin oscillations [8], whereas bumps in the 3 MeV 
region of rare earth nuclei are now determined to be of 
Ml character [9,10]. The electromagnetic character and 
measured strength of the latter pygmy resonance is com- 



patible with the scissors mode. Recently [11,12], the RSF 
picture of iron isotopes has been further modified by the 
observation of an anomalous increase in strength at 7 en- 
ergies below 4 MeV. It is clear that the present situation 
urges for new experimental results. 

The stable molybdenum isotopes are well suited as 
targets for the study of nuclear properties when going 
from spherical to deformed shapes. In this work we 
perform a systematic analysis of the RSFs of the six 
93_98 Mo isotopes. The RSFs depend on the dynamic 
properties of electric charges present within these sys- 
tems (Z — 42). Since the nuclear deformation varies 
from spherical shapes [fi ~ 0) at iV = 51 to deformed 
shapes (j3 ~ 0.2) at N = 56, we expect to observe effects 
due to shape changes. Furthermore, these nuclei reveal 
weak GEDR tails at low E 1 making them interesting ob- 
jects in the search for other weak structures in the RSF. 

The Oslo Cyclotron group has developed a sensitive 
tool to investigate RSFs for E 1 below the neutron binding 
energy S n . The method is based on the extraction of pri- 
mary 7-ray spectra at various initial excitation energies 
Ei measured in particle reactions with one and only one 
charged ejectile. From such a set of primary 7-spectra, 
nuclear level densities and RSFs can be extracted [13-15]. 
The level density reveals essential nuclear structure infor- 
mation such as thermodynamic properties and pair cor- 
relations as function of temperature. These aspects of 
the molybdenum isotopes will be presented in a forth- 
coming work. Various applications of the Oslo method 
have been described in Refs. [16-23]. 
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II. EXPERIMENTAL METHOD 

The particle-7 coincidence experiments were carried 
out at the Oslo Cyclotron Laboratory for 93_98 Mo using 
the CACTUS multi-detector array. The charged ejectiles 
were detected with eight particle telescopes placed at an 
angle of 45° relative to the beam direction. An array of 
28 Nal 7-ray detectors with a total efficiency of <~15% 
surrounded the target and particle detectors. 

In the present work, results from eight different re- 
actions on four different targets are discussed. Results 
from two of those reactions have been reported earlier. 
The beam energies for the different reactions are given 
in parentheses: 

1) 94 Mo( 3 Hc,a 7 ) 93 Mo (new, 30 McV) 

2) 94 Mo( 3 He, 3 He'7) 94 Mo (new, 30 MeV) 

3) 96 Mo( 3 He,a7) 95 Mo (new, 30 McV) 

4) 96 Mo( 3 He, 3 He' 7 ) 96 Mo (new, 30 MeV) 

5) 97 Mo( 3 He,crr) 96 Mo (reported in [11,23], 45 MeV) 

6) 97 Mo( 3 He, 3 He'7) 97 Mo (reported in [11,23], 45 MeV) 

7) 98 Mo( 3 He,a7) 97 Mo (new, 45 MeV) 

8) 98 Mo( 3 He, 3 He'7) 98 Mo (new, 45 MeV) 

The targets were self-supporting metal foils enriched to 
~ 95% with thicknesses of ~ 2 mg/cm 2 . The experiments 
were run with beam currents of ~ 2 nA for 1-2 weeks. 
The reaction spin windows are typically 1^2 — 6h. 

The experimental extraction procedure and the as- 
sumptions made are described in Refs. [13,15] and refer- 
ences therein. For each initial excitation energy Ei, de- 
termined from the ejectile energy and reaction Q value, 
7-ray spectra are recorded. Then the spectra are un- 
folded using the known 7-ray response function of the 
CACTUS array [24] . These unfolded spectra are the ba- 
sis for making the first generation (or primary) 7-ray ma- 
trix [25] , which is factorized according to the Brink- Axel 
hypothesis [4,26] as 



P(Ei, £" 7 ) oc p(Ei — E 1 )T{E 1 



(1) 



Here, p is the level density and T is the radiative trans- 
mission coefficient. 

The p and T functions can be determined by an itera- 
tive procedure [15] through the adjustment of each data 
point of these two functions until a global x 2 minimum of 
the fit to the experimental P{Ei,E 1 ) matrix is reached. 
It has been shown [15] that if one solution for the multi- 
plicative functions p and T is known, one may construct 
an infinite number of other functions, which give identical 
fits to the P matrix by 

p{E { - £ 7 ) = Acxp[a(E t - E 7 )} p(E t - £ 7 ), (2) 
f(E 7 ) = Bexp(aE y )T(E y ). (3) 

Consequently, neither the slope (a) nor the absolute val- 
ues of the two functions (A and B) can be obtained 
through the fitting procedure. 



The parameters A and a can be determined by nor- 
malizing the level density to the number of known dis- 
crete levels at low excitation energy [27] and to the level 
density estimated from neutron-resonance spacing data 
at the neutron binding energy S n [28]. The procedure 
for extracting the total level density p from the reso- 
nance energy spacing D is described in Ref. [15]. Here, 
we will only discuss the determination of parameter B of 
Eq. (3), which gives the absolute normalization of T. For 
this purpose we utilize experimental data on the average 
total radiative width of neutron resonances at S n (T 7 ). 

We assume here that the 7-decay in the continuum is 
dominated by El and Ml transitions. For initial spin / 
and parity tt at S n , the width can be written in terms of 
the transmission coefficient by [29] 



(T T > 



p(S n -E 7 ,I f ), (4) 



where the summation and integration run over all final 
levels with spin //, which are accessible by 7 radiation 
with energy E 1 and multipolarity El or Ml. 

A few considerations have to be made before B can 
be determined. Methodical difficulties in the primary 7- 
ray extraction prevents determination of the functions 
T(i? 7 ) in the interval E 7 < 1 MeV and p(E) in the 
interval E > S n - 1 MeV. In addition, T(E 1 ) at the 
highest 7-energies, above E 7 ~ S n — 1 MeV, suffers from 
poor statistics. For the extrapolation of p we apply the 
back-shifted Fermi gas level density as demonstrated in 
Ref. [22] . For the extrapolations of T we use an exponen- 
tial form. As a typical example, the extrapolations for 
98 Mo are shown in Fig. 1. The contribution of the ex- 
trapolations of p and T to the calculated radiative width 
in Eq. (4) does not exceed 15% [20]. The experimental 
widths (r 7 ) in Eq. (4) arc listed in Table I. For 94 Mo, this 
width is unknown and is estimated by an extrapolation 
based on the 96 Mo and 98 Mo values. 

The total radiative strength function for dipole radi- 
ation (L = 1) can be calculated from the normalized 
transmission coefficient T by 



f(E 7 ) 



1 T(£7 7 ) 



2tt 



i? 3 



(5) 



The RSFs extracted from the eight reactions are dis- 
played in Fig. 2. As expected, the RSFs seem not to 
show any clear odd-even mass differences. The results 
obtained for the ( 3 He,a) and ( 3 He, 3 He') reactions popu- 
lating the same residual nucleus reveal very similar RSFs. 
Also for 96 Mo two different beam energies have been ap- 
plied, giving very similar RSFs. Thus, the observed en- 
ergy and reaction independency gives further confidence 
in the Oslo method. 
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III. DESCRIPTION OF THE RADIATIVE 
STRENGTH FUNCTIONS 



An inspection of the experimental RSFs of Fig. 2 re- 
veals that the RSFs are increasing functions of 7 energy 
for E n > 3 MeV. This indicates that the RSFs are in- 
fluenced by the tails of the giant resonances. As follows 
from previous work, the main contribution (about 80%) 
is due to the electric dipole resonance (GEDR). The mag- 
netic resonance (GMDR) and the isoscalar E2 resonance 
are also present in this region. 

If the GEDR is described by a Lorentzian function, 
one will find that the strength function approaches zero 
in the limit E 1 — > 0. However, the 144 Nd(n,7a) reac- 
tion [31] strongly suggests that /ei has a finit value in 
this limit. Kadmenskh, Markushev and Furman (KMF) 
have developed a model [32] describing this feature for 
the electric dipole RSF: 



_l o.7 < 7 E1 r| 1 (^ + 47r 2 r 2 ) 

3^ 2 c 2 £ E i(£2 - £| i)2 



(6) 



The temperature T depends on the final state / and for 
simplicity we adapt the schematic form 



T(E f ) = JUfJa, 



(7) 



where the level density parameter is parametrized as 
a = 0.21^4°' 87 MeV -1 . The intrinsic energy is estimated 
by Uf = Ef — C\ — -Epair with a back-shift parameter 
of" Ci = -6.6A~ - 32 MeV [33]. The pairing energy con- 
tribution impair is evaluated from the three point mass 
formula of Ref . [35] . 

Although this model has been developed for spheri- 
cal nuclei, it has been successfully applied to 56 ' 57 Fe and 
several rare earth nuclei [12,20-22] assuming a constant 
temperature parameter T in Eq. (6), i.e. one that is in- 
dependent of excitation energy. In this work we assume 
that the temperature depends on excitation energy ac- 
cording to Eq. (7), which gives an increase in the RSF at 
low 7 energy. 

In order to compare the excitation energy dependent 
RSF with experiments, Eq. (6) should be folded with the 
final excitation energies of the specific experiment giving 



</ei(£ 7 )> 



E, 



±— \^ El dE f f E1 (E f ), (8) 



where the integration runs over all final excitation en- 
ergies Ef which are experimentally accessible by tran- 
sition with given energy E 1 . Our data set of primary 
7-ray spectra is typically taken between E\ ~ 4 MeV 
and Ei ~ S n initial excitation energy. 

The Ml GMDR contribution to the total RSF is de- 
scribed by a Lorentzian form through 



/M1 ^ j ~3^V {Ei-E^Y+E^ 



(9) 



This approach is in accordance with numerous experi- 
mental data obtained so far, and is recommended in Ref. 
[28]. 

Although of minor importance, the E2 radiation 
strength /e2 has also been included. Here, we use the 
same Lorentzian function as in Eq. (9), but with different 
resonance parameters and an additional factor 3/(5i? 2 ). 
The resonance parameters for the El, Ml and E2 res- 
onances are taken from the compilation of Refs. [28,34] 
and are listed in Table I. 

The enhanced RSF at low 7 energies has at present 
no theoretical explanation. Recently, the same enhance- 
ment has been observed in the iron isotopes [11,12]. We 
call this structure a soft pole in the RSF and choose a 
simple power law parametrization given by 



1 



softpolc 



3rfc 2 



AE~ 



(10) 



where A and b are fit parameters, and E 1 is given in 
MeV. 

Previously, a pygmy resonance around E 7 ~ 3 MeV 
has been reported in several rare earth nuclei [20-22]. 
The electromagnetic character of the corresponding RSF 
structure is now established to be of Ml type [9,10] and 
is interpreted as the scissors mode. Deformed nuclei can 
in principle possess this collective motion, and e.g., 98 Mo 
with a deformation of (i ~ 0.18, could eventually show 
some reminiscence of the scissors mode. Data on 94 Mo 
[36] and 96 Mo [37] show a summed Ml strength to mixed 
symmetry 1 + states around ^3.2 MeV in the order of 
~ 0.6/^. This is about one order of magnitude lower 
than the Ml strength observed in well-deformed rare 
earth nuclei using the present method. This Ml strength 
is deemed too weak to cause a visible bump in our RSFs 
above 3 MeV. 

We conclude that a reasonable composition of the total 
RSF is 



/ — k(/e1 + /Ml + /softpolc) + i^/E2, 



(11) 



where k is a normalization constant. Generally, its value 
deviates from unity for several reasons; the most impor- 
tant reasons are theoretical uncertainties in the KMF 
model and the evaluation of B in Eq. (4). We use k, 
A and b as free parameters in the fitting procedure, and 
the results for the eight reactions are summarized in Ta- 
ble II. 

In Fig. 3 the various contributions to the total RSF of 
98 Mo are shown. The main components are the GEDR 
resonance and the unknown low energy structure. We ob- 
serve that the El component exhibits an increased yield 
for the lowest 7 energies due to the increase in temper- 
ature T. However, this effect is not strong enough to 
explain the low energy upbend. 

Figure 2 shows the fit functions for all reactions and 
gives qualitative good agreements with the experimental 
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data. The fitting parameters k, A and b are all sim- 
ilar within the uncertainties. It should be noted that 
the soft pole parameters coincide with the description of 
the 57 Fe nucleus [12] having A = 0.47(7) mb/MeV, and 
6 = 2.3(2). 

The almost a factor-of-two higher values of the RSFs 
for E 1 > 3 MeV when going from TV = 51 to 56 can 
be understood from the corresponding evolution of nu- 
clear deformation. Following the onset of deformation 
the GEDR will split into two parts, where 1/3 of its 
strength is shifted down in energy and 2/3 up. Pho- 
toneutron cross sections [34] do not show a splitting into 
two separate bumps, however, the observed increase in 
width Tei as function of neutron number (see Table I) 
supports the idea of a splitting, which is well known from 
other more deformed nuclei. Figure 2 demonstrates that 
the adopted widths describe very well the variation of the 
RSF strength as function of mass number. 

In order to test that the prominent soft pole struc- 
ture is present in the whole excitation energy region, we 
have performed the following test. Assuming that the 
level density from Eq. (1) is correct, we can estimate the 
shape of the strength functions starting at various initial 
excitation energies using 



f(E 7 ,Ei) — 



1 JV(E i )P(E l , E 7 ) 
2tt p(Ei-EJE* 



(12) 



Actually, f{E 1 ,Ej) would have been the proper expres- 
sion to investigate, but due to technical reasons we choose 
f{E~ l ,Ei). One problem is that the normalization con- 
stant is only roughly known through the estimate 



M{Ei) 



jf dE lP (E t - E,)T(E 7 ) 
Jo % dE 7 P(Ei, E y ) 



(13) 



with Ei < S n - However, for the expression /(i? 7 , £^) we 
are only interested in the shape of the RSFs, and an ex- 
act normalization is therefore not crucial. The evaluation 
assumes that eventual temperature dependent behavior 
of the RSF is small compared to the soft pole structure 1 . 

In Fig. 4, the RSFs for 96 < 98 Mo are shown at various 
initial energies Ei. For comparison, the figure also in- 
cludes the global RSFs (solid lines) obtained with the 
Oslo method (Fig. 2). Within the error bars the data 
support that the soft pole is present in all the excitation 
bins studied. 

The origin of the soft pole cannot be explained by any 
known theoretical model. One would therefore need to 
know the 7-ray multipolarity as guidence for theoretical 
approaches to this phenomenon. Rough estimates of the 
reduced strength can be obtained from 



B(XL T) 



1 L[(2L+ l)!!] 2 / he 



8tt 3 



I. 



(14) 



The integrated cross-section / is determined by 

r 3MoV 

d£ 7 o-(£; 7 ), 

' lMcV 

and the average energy is given by 

r 3MeV 



Jib 



(15) 



where 



(E^^r 1 dE 7 E 7 a(E 7 ), 

JlMoV 



a(E 7 ) = 3(nhc) 2 E 1 f so { tpo i e (E 1 ). 



(16) 



(17) 



In the evaluation, we have integrated the soft pole be- 
tween 1 and 3 MeV. Thus, the estimates give a lower limit 
for the respective B(XL j) values. The correct result will 
of course depend on the functional form of /softpok^-E^y) 
below 1 MeV, however, no experimental data exist in this 
region and any assumption here would be highly specu- 
lative. 

With the assumptions above, we get on the average a 
B{El) value of 0.02 e 2 fm 2 , which is 0.07 % of the sum 
rule for the GEDR. For B(M1), we get roughly 1.8(1%, 
which is 3 — 4 times larger than the observed strength 
to mixed symmetry 1 + states around 3 MeV [36,37]. Fi- 
nally, we obtain a B(E2) value around 8000c 2 fm 4 , which 
is 3 — 8 times larger that for the excitation to the first ex- 
cited 2 + states in the even molybdenum isotopes. Thus, 
we cannot exclude any of these multipolarities, since nei- 
ther of them would yield unreasonably high transition 
strengths. 



IV. SUMMARY AND CONCLUSIONS 



As expected, the observed RSFs reveal very similar 
shapes since they all refer to isotopes with the same nu- 
clear charge. When going from N = 51 to 56 the RSF 
increases by almost a factor of two for E 7 > 3 MeV, 
which can be understood from the change of nuclear de- 
formation. With the onset of deformation, the increasing 
resonance GEDR width Tei is responsible for the increas- 
ing strength. 

An enhanced strength at low 7 energies is observed, 
which is equally strong for all isotopes and excitation 
energies studied. A similar enhancement has also been 
seen in the iron isotopes. The multipolarity of the soft 
pole radiation is unknown and there is still no theoretical 
explanation for this very interesting phenomenon. 



1 Simulations using the KMF model indicate a maximum 20 % effect from temperature dependence of the RSF. 
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TABLE I. Parameters used for the radiative strength functions. The data are taken from Ref. [28]. The El resonance 
parameters for the even Mo isotopes are based on photo absorption experiments [34], and the parameters for the odd Mo 
isotopes are derived from interpolations. 



Nucleus 


Eei 


C"E1 


Tei 


-Emi 


<TM1 




E E 2 


CTE2 


Te2 


<r 7 ) 
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93 Mo 


16.59 


173.5 


4.82 


9.05 


0.86 


4.0 


13.91 


2.26 


4.99 


160(20) 


94 Mo 


16.36 


185.0 


5.50 


9.02 


1.26 


4.0 


13.86 


2.24 


4.98 
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4.0 
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97 Mo 


16.00 


187.0 


5.98 


8.92 


1.58 


4.0 


13.71 


2.19 


4.95 


110(15) 


98 Mo 


15.80 


189.0 


5.94 


8.89 


1.65 


4.0 


13.66 


2.17 


4.93 


130(20) 



"Estimated from systematics. 



TABLE II. Parameters obtained from the fitting of the radiative strength function. 



Reaction 


A 

(mb/MeV) 
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( 3 He,a) 93 Mo 


0.37(7) 


2.6(3) 


0.44(4) 


( 3 He, 3 He') 94 Mo 


0.48(5) 


2.5(2) 


0.36(2) 


( 3 He,a) 95 Mo 


0.48(6) 


2.6(2) 


0.39(2) 


( 3 He, 3 He') 96 Mo 


0.60(4) 


3.2(2) 


0.36(1) 


( 3 He,a) 96 Mo 


0.47(14) 


2.7(6) 


0.32(4) 


( 3 He, 3 He') 97 Mo 


0.47(7) 


2.4(3) 


0.38(3) 


( 3 He,a) 97 Mo 


0.30(10) 


2.2(5) 


0.45(5) 


( 3 He, 3 He') 98 Mo 


0.22(7) 


2.1(5) 


0.52(4) 
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FIG. 1. Measured level density p (upper panel) and radiative transmission coefficient T (lower panel) for Mo. The straight 
lines are extrapolations needed to calculate the normalization integral of Eq. (4) . The triangle in the upper panel is based on 
resonance spacing data at S n - 
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FIG. 2. Normalized RSFs for 93 98 Mo. The filled and open circles represent data taken with the ( 3 He,a) and ( 3 He, 3 He') 
reactions, respectively. The filled triangles in 93 ' 95 Mo are estimates of El RSF of hard primary 7-rays [30] . The solid and 
dashed lines are fits to the RSF data from the two respective reactions (see text). 
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FIG. 3. Experimental radiative strength function of Mo compared to a model description including GEDR, GMDR and 
the isoscalar E2 resonance. The empirical soft pole component is used to describe the low energy part of the RSF. 
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FIG. 4. RSFs for 96 ' 98 Mo at various initial excitation energies. The soft pole is present for all Ei. The solid lines display the 
RSFs obtained in Fig. 2. 
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